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Abstract 



Measurements of event-by-event fluctuations of the mean transverse momentum in 
Pb-Au collisions at 40, 80, and 158 A GeV/c are presented. A significant excess of 
mean px fluctuations at mid-rapidity is observed over the expectation from statisti- 
cally independent particle emission. The results are somewhat smaller than recent 
measurements at RHIC. A possible non-monotonic behaviour of the mean pj- fluc- 
tuations as function of collision energy, which may have indicated that the system 
has passed the critical point of the QCD phase diagram in the range of (is under 
investigation, has not been observed. The centrality dependence of mean pt fluctu- 
ations in Pb-Au is consistent with an extrapolation from pp collisions assuming that 
the non-statistical fluctuations scale with multiplicity. The results are compared to 
calculations by the RQMD and URQMD event generators. 

Key words: Nuclear reactions 197 Au(Pb, X), E = 40, 80, 158,4 GeV; 
event-by-event transverse momentum fluctuations, QCD phase transition, critical 
point. 

PACS: 25.75.Gz, 25.75.Nq 



1 Introduction 

The investigation of high energy nucleus-nucleus collisions provides a unique 
tool to study the properties of hot and dense nuclear matter. The motivation 
largely comes from QCD calculations, particularly on the lattice, which pre- 
dict, at sufficiently high temperatures and densities, a transition from hadronic 
matter to a plasma of deconfined quarks and gluons. In this phase chiral sym- 
metry is also restored. For the case of vanishing net baryon density (or baryon 
chemical potential fin), state of the art lattice calculations suggest that the 
QCD phase transition occurs at a critical temperature of T c =175±8 MeV [1]. 
However, a full exploration of the QCD phase diagram in the T-/i B -plane 
is desired to obtain insight into the mechanism of deconfinement and chi- 
ral symmetry restoration, as well as the properties of astrophysical objects 
such as neutron stars. At finite net baryon chemical potential hb, lattice cal- 
culations run into technical limitations, which have been partially overcome 
recently [2,3], showing how T c drops with increasing /j, B . Complex structures 
of the QCD phase diagram have been unveiled, such as the possible existence 
of a critical point of the phase boundary, and a color-superconducting phase 
at very high baryon density and low temperature [2,4,5,6]. 
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Most theoretical studies of the phase transition imply that it is second or- 
der or a continuous but rapid cross-over, at least for values less than 
500 MeV [1,2,3]. Passage of a system through a second order transition or 
close to a critical point may lead to critical phenomena, long-range correla- 
tions and large fluctuations. The study of event-by-event fluctuations therefore 
provides a novel probe to explore the QCD phase diagram, searching for the 
quark-gluon plasma (QGP) and the QCD critical point. Such measurements 
became possible with large acceptance experiments at SPS and RHIC, where 
the high multiplicity of charged particles produced in collisions of lead and 
gold nuclei allows a precise determination of global observables on an event- 
by-event basis. Pioneering experimental studies have been performed by the 
NA49 experiment at the SPS, on fluctuations of the mean transverse momen- 
tum p T [7] and of the K/n ratio [8] in central Pb-Pb collisions at 158 A GeV/c. 

Mean p T fluctuations have been investigated in pp collisions at the ISR [9]. 
These data may be used as a reference to study the degree of thermal equili- 
bration in nucleus- nucleus collisions [10,11]. On the other hand, an enhance- 
ment of mean p? fluctuations might be related to non-trivial effects showing 
up specifically in nuclear reactions. It was predicted that mean p T fluctua- 
tions can be enhanced if the system passes through the QCD critical point, 
where long wave length fluctuations of the sigma field develop, leading to 
fluctuations of pions through the strong cr-n-Tr coupling [12,13]. Decay of the 
Polyakov loop [14] and random color fluctuations [15] are also proposed as pos- 
sible mechanisms to enhance momentum fluctuations when the system passes 
the phase boundary from the QGP to the hadronic phase. At high baryon 
chemical potential ~ 3T), the possible occurence of large baryon number 
density fluctuations has been discussed [16]. 

The NA49 experiment observed non-statistical mean px fluctuations consistent 
with zero [7] at forward rapidity in 158 A GeV/c Pb-Pb collisions at the SPS. 
This result was explained by a small positive contribution from Bose-Einstein 
correlations, which is quantitatively compensated by a negative contribution 
arising from the finite two-track resolution of the detector. No indications 
for additional non-statistical fluctuations were observed at forward rapidity. 
However, the complex structure of the QCD phase diagram requires a detailed 
study of event-by-event fluctuations in the central rapidity region and at all 
available beam energies. In particular, the data taken at lower SPS energies 
allow for an investigation of event-by-event fluctuations at larger The 
observation of a possible non-monotonic behaviour of the fluctuation strength 
as function of \1b may indicate the location of the critical point in the QCD 
phase diagram. 

In this paper, we present measurements of event-by-event fluctuations of the 
mean transverse momentum px of charged particles produced near mid-rapidity 
in Pb-Au collisions at 40, 80, and 158 A GeV/c. 
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Fig. 1. The CERES experiment. 
2 Experimental setup 

The CERES spectrometer at the CERN-SPS was originally conceived for the 
measurement of low mass e + e~-pairs produced in proton and heavy-ion in- 
duced collisions with nuclei [17]. This version of CERES consisted of two ra- 
dial silicon drift detectors (SDD's) and two ring-imaging Cherenkov counters 
(RICH1,2), supplemented with a pad detector behind the RICHes. In 1998, 
after removal of the pad detector, a major upgrade of the spectrometer was 
performed by the addition of a large cylindrical Time Projection Chamber 
(TPC) (see Fig. 1). All subdetectors have a common acceptance in the polar 
angle range of 8° < 6 < 15° at full azimuth, corresponding to a pseudorapidity 
acceptance of 2.1 < rj < 2.65. 

The SDD's are located about 12 cm downstream of the target system. Each 
SDD has a uniform radial drift field, with 360 equally spaced readout anodes 
arranged along the outer perimeter of the Si-wafer. A SDD track segment is 
defined by a hit in each of the two detectors whose straight line connection 
points to the vertex of the interaction. 

The TPC, located about 3.8 m downstream of the target, has an active length 
of 2 m and a diameter of 2.6 m. It is filled with a mixture of 80% Neon and 
20% C0 2 , and operated inside the magnetic field generated by two opposite- 
polarity solenoidal coils, which are placed around the TPC. In the TPC drift 
volume, the ionization electrons created along charged particle tracks drift 
outwards along the field lines of a radial electric drift field. They are detected 
in 16 multi-wire proportional chambers with cathode pad readout, which form 
the polygon- shaped outer termination of the active TPC volume. The TPC 
provides the reconstruction of up to 20 space points along a curved charged 
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particle track, thereby improving the momentum resolution of the spectrom- 
eter and supplementing its particle identification capability via the measure- 
ment of the specific energy loss dE/dx. 

In 1999 (40 A GeV/c Pb-beam period) the target system consisted of eight Au 
foils of 25 /mi thickness, separated by 3.1 mm in beam direction and adding 
to 0.82% of a hadronic interaction length. In 2000 (80 and 158 A GeV/c Pb- 
beam period) the target was replaced by 13 Au foils of the same dimensions 
and 2 mm spacing, with a total thickness of 1.33% of an interaction length. 

The beam trigger (BEAM) is defined by the coincidence of two Beam Counters 
BC1 and BC2 located in front of the target. Absence of a signal in a third Beam 
Counter BC3 located behind the target, in coincidence with BEAM, defines 
the interaction trigger (INT). An online centrality selection was performed by 
an additional threshold on the integrated pulse height in the SDD in 1999, or 
on the signal in a scintillation Multiplicity Counter (MC) in 2000. The data 
taking rate was 300-500 events per burst, at a typical Pb beam intensity of 
10 6 ions per burst. In 1999, the TPC readout was not yet fully operational. As 
a consequence, the signals of only 60% of the 15360 TPC channels could be 
stored. The 80 A GeV/c data were taken while setting up for the 158 A GeV/c 
beam period, therefore no SDD information is available for the 80 A GeV/c 
data. 



3 Data analysis 

3.1 Event selection and track reconstruction 

In the present analysis we have applied an offline centrality selection of the 
upper 20% of the total geometric cross section. The centrality of the collision 
was defined by the number of charged particle tracks reconstructed in the SDD 
(for the 40 A GeV/c data) or the pulse height in the Multiplicity Counter 
MC (for the 80 and 158 A GeV/c data). The corresponding mean number of 
participating nucleons (N pait ) and mean number of nucleon-nucleon collisions 

(N co \\) was derived from a geometric nuclear overlap model using b = \Jo~fn 
and resulting, with gnat = 30 mb, in a total cross section of o~g — 6-94 b [18]. 
Our classification of central events comprises the 6.5% most central fraction 
of the total geometric cross section. For the centrality dependent studies we 
have subdivided our sample into the 0-5%, 5-10%, 10-15%, and 15-20% most 
central events (see Table 1). 

Additional offline cuts have been applied to reject pile-up events and non- 
target interactions. The total number of analyzed events, after centrality and 
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Table 1 

Definition of centrality classes. 



a/(Tg eo 


- 6.5% 


0-5% 


5 - 10% 


10 - 15% 


15 - 20% 




337 


347 


287 


239 


200 


(iVcoll) 


781 


809 


638 


503 


393 



quality cuts, is 1.5M at 40 A GeV/c and 0.5M at 80 and 158 A GeV/c, 
respectively. 

SDD track information has only been used for the centrality definition in 
40 A GeV/c events, and for a number of systematic consistency checks as 
described in Section 3.3. All final results presented below are based on TPC 
track information only. 

A TPC track is defined by a number of hits reconstructed in subsequent 
readout-planes of the TPC and matching a momentum-dependent track model 
generated by a Monte Carlo simulation. Depending on the polar angle, a TPC 
track consists of up to 20 hits. The momentum of the track is evaluated from 
the azimuthal deflection of the track inside the active TPC volume. At the 
stage of the calibration used for this analysis, we achieved a momentum reso- 
lution of about Ap/p = ((1.5%) 2 + (1.4%-p (GeV/c)) 2 ) 1 / 2 for tracks containing 
more than 18 hits (in the 80 and 158 A GeV/c data sets). 

A number of fiducial and quality cuts have been applied to provide stable 
tracking conditions and to reject tracks from secondary particles: 

• In the 40 A GeV/c data set, a fiducial cut in the azimuthal angle 17/247T < 
< 27r is applied. Outside this region the efficiency was low due to mal- 
functioning read-out electronics during the 1999 beam period. No such cut 
was applied for the 80 and 158 A GeV/c data from the year 2000, where 
the read-out worked properly. The effect of the azimuthal fiducial cut on 
the final results has been studied and is included in the systematic errors 
(see Section 3.3). 

• For the analysis of the mean px fluctuations only TPC tracks inside the 
full-length TPC track acceptance (2.2 < rj < 2.7) were used to provide 
sufficient px resolution. For the same reason, the analysis was restricted to 
tracks with transverse momenta 0.1 < p T < 1-5 GeV/c. 

• The minimum number of hits per track is 11 at 40 A GeV/c and 12 at 80 
and 158 A GeV/c in the full-length TPC track acceptance 2.2 < rj < 2.7. 

• To suppress secondary particles it is required that the back-extrapolation 
of the particle trajectory into the target plane misses the interaction point 
by no more than 10 cm in transverse direction. 
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• The x 2 I dof of the track momentum fit has to be less than three times the 
r.m.s. of the x 2 1 dof distribution. 

These cuts are used for all final results presented below. A systematic variation 
of these cuts has been applied to estimate our systematic uncertainties, as 
described in Section 3.3. 



3.2 Measures of mean px fluctuations 



In the following, we briefly review quantities which have been proposed as 
measures for event-wise mean pr fluctuations and summarize notations used 
in this paper. 

The event-wise mean M x of a single-particle observable x, averaged over par- 
ticles i in the acceptance of an event j, is given by: 

Mi = (1) 



where Nj is the multiplicity of event j. 

For a quantity Xj which is defined for each event j we calculate its mean over 
all events: 

(X) = E ^ lV0 , (2) 
n 



where n indicates the total number of events. The weighting factor Wj for each 
event j is equal to Nj/ (N) for variables Xj which are already an average over 
Nj particles, i.e. Xj = M 3 X . For other quantities Xj, such as multiplicity Nj, 
the weighting factor Wj is equal to 1. 

The variance of the distribution of X is given by: 

(AX 2 ) = (X 2 ) - (X) 2 . (3) 



The inclusive mean (the mean over all particles in all events) and the variance 
of a single particle observable x are defined as: 
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and 



Ax 2 = x 2 ~x 2 . (5) 

The mean and the variance of M x are obtained by substituting M x for X in the 
above equations, and including the event multiplicities as appropriate weight- 
ing factors. This weighting procedure provides the most precise estimate of the 
variance of the parent distribution in case of a finite mean multiplicity [13]: 

„ =* (6) 



and 



In the present event-by-event analysis, we search for dynamical mean p T fluc- 
tuations beyond those expected in a purely statistical scenario. Dynamical 
mean p? fluctuations would therefore result in an event-by-event distribution 
of M PT which is wider than that expected from the inclusive pt distribution 
and the finite event multiplicity. 

In previous analyses, the measure $ PT has been used to quantify non-statistical 
mean p T fluctuations [10]: 



$ = 



{Z2) (8) 



\| (N) 



N 

where z and Z are defined as z = px — Pt for each particle, and Z = Zj for 

i=i 

each event, respectively. The measure <3> PT vanishes in the absence of correla- 
tions and dynamical mean p? fluctuations. There is an approximate expression 
for $ PT in terms of the variances of the event- wise M PT and the r.m.s. of the 
inclusive Pt distributions [19]: 



% T = \/( AM pt) ~ V^4- (9) 

A different measure for dynamical mean p T fluctuations has been proposed 
in [19]: 
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This expression provides a direct relation between the variance of the inclusive 
Pt distribution, the mean multiplicity, and the variance of the event-by-event 
mean p T distribution. In case of vanishing non-statistical fluctuations and 
correlations, <7 2 VT ^ yn is equal to zero. It has also been shown that o 2 VT ^ yn is the 
mean of the covariances of all possible pairs of two particles from the same 
event [19]. 

An alternative approach to evaluate o^ Tdyn is based on the analysis of sub- 
events [19]. The subevent method was applied for a number of consistency 
checks, however, all final results presented below are obtained using Eq. (10) 
for the evaluation of o^ Tdyn . 

There is an important relation between $ PT and ai T d [19]: 



.2 2$ PT y/Ap| 

p T ,dyn i A t\ 



Based on our results presented below we find that this relation holds extremely 
well. 

In order to account for a possible change of mean pr at different beam ener- 
gies, we define a dimensionless measure, the "normalized dynamical fluctua- 
tion" XL ■ 



v — (2 \ V a PT,dyn , 1oA 

S PT = sgn{a pT4yn ) ■ = . (12) 



Due to the central limit theorem, the two measures S PT and $ PT are equal to 
zero in a purely statistical distribution. If dynamical mean pt fluctuations are 
present, both measures are finite and positive. However, also (anti-) correla- 
tions in momentum space cause non- vanishing values. Long-range correlations 
occur as a consequence of energy and momentum conservation, while Bose 
and Fermi statistics, final state interactions, and experimental effects such as 
the finite two-track resolution are the origin of short-range correlations. Since 
such correlations mask the true fluctuation signal, it is crucial to investigate 
their contribution to Eq. (8) and Eq. (10) quantitatively, as will be described 
in Section 4. 

The best way for a quantitative determination of mean p T fluctuations is 
clearly the use of a dimensionless measure. In addition, a comparison between 
different experiments should be possible. In this context we note that <3> PT is 
neither dimensionless nor independent of the event multiplicity. This makes a 
comparison between experiments and to theory difficult because the multiplic- 
ity depends on the acceptance window of the experiment and on beam energy. 
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Fig. 2. The mean charged particle multiplicity (N) (left) and the fluctuation 
measures & PT (middle) and E PT (right) in central Pb-Au events at 40, 80, and 
158 A GeV/c as function of the pseudorapidity bin size Ar/. The center of the At? 
window is always fixed at r\ = 2.45. Note that the statistical errors are correlated. 

Since different contributions to the fluctuation signal have different multiplic- 
ity dependences, the multiplicity dependence of $ PT is a priori unknown. 

In Fig. 2 (left panel) the mean multiplicity is shown as function of the size of 
the pseudorapidity acceptance window Ar/ at 40, 80, and 158 A GeV/c. The 
much smaller multiplicity at 40 A GeV/c is related to the smaller azimuthal 
fiducial window in this data set. The measure $ PT , shown in Fig. 2 (middle 
panel), exhibits a strong dependence on Ar/ and beam energy. However, com- 
parison to the left panel in Fig. 2 indicates that this dependence is strongly 
correlated with the change in mean multiplicity. In contrast, the normalized 
dynamical fluctuation S PT changes very little with beam energy, as shown in 
Fig. 2 (right panel). The increase of S PT towards small Ar/ may be attributed 
to physical correlations occuring on small Ar/ scales. These contributions ap- 
pear to vanish at Ar/ > 0.5 and S PT does not any longer depend on Ar/, 
although the multiplicity is still increasing 1 . This suggests that S PT is a more 
universal measure of fluctuations and it will be mainly used in the present 
work. In Section 4 we will discuss the implications of the results presented in 
this figure. 



1 We have checked that the definition of at ,,„„ in terms of the mean of covariances 

p T ,dyn 

and the definition in Eq. (10) lead to the same results for S PT within statistical 
errors, also for small multiplicities. In particular, the dependence on Ar/ is observed 
using both definitions. 
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3.3 Study of systematic uncertainties 



In this subsection, we briefly summarize the different contributions to the sys- 
tematic uncertainties of the fluctuation measurements presented below. The 
contributions to the systematic errors of the mean pt fluctuation measure- 
ments are listed in Table 2. 

Biases due to our analysis chain have been evaluated with the help of the 
URQMD (version 1.1) event generator [20] (without rescattering) , which ex- 
hibits mean fluctuations of similar magnitude as observed in the data. The 
generated events were processed through a full Monte Carlo simulation of the 
detector using the detector simulation package geant [21] including signal 
digitization, and analyzed by our reconstruction software applying all cuts 
listed before. No significant bias on the initial fluctuation signal was observed 
after event reconstruction. This implies also that a small remaining contribu- 
tion of secondary particles does not alter the fluctuation signal. 
The systematic uncertainty due to the finite tracking efficiency was evaluated 
using the same Monte Carlo chain as before. The number reported in Table 2 
corresponds to the change of the reconstructed S PT when randomly a fraction 
of the tracks (up to 20%) is artificially removed from the sample. The esti- 
mated tracking efficiency in real data is 90%. 

An online before- and after-protection of ±1 /is against a second beam parti- 
cle was applied during data taking. In the offline analysis, we systematically 
increased this protection to ±3 /is via event selection but found very little 
effect on the reconstructed value of S PT . 

The absolute momentum scale of the spectrometer has been verified by the 
peak position of the invariant mass distributions of reconstructed A and K° 
decays [22]. At the present stage of the calibration, the uncertainty of the 
absolute momentum scale is at most ±0.02 (GeV/c) -1 in l/p. We have arti- 
ficially shifted l/p of each reconstructed track by a constant offset to study 
the effect of this uncertainty. As a test of the long-term stability of the mo- 
mentum scale we divided our events into subsamples which were taken close 
in time and analyzed them separately. After proper calibration the results of 
the subsamples were consistent with each other and with the result of the full 
data set. 

The impact of the azimuthal fiducial cut in the 40 A GeV/c data set was 
studied by applying the same cut to the 158 A GeV/c data and comparison 
to the full acceptance result. We found that the contribution of the fiducial 
cut to the final result is very small. 

The inclusion of SDD track information leads to a powerful rejection of non- 
vertex tracks, if only TPC tracks with a match to the SDD are used in the 
analysis. We compared the fluctuation results with and without use of the 
SDD in the 158 A GeV/c data set. From the difference we estimated the sys- 
tematic uncertainty if SDD information is not used. Since SDD information 
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Table 2 

Systematic errors of S PT at 2.2 < 77 < 2.7 in the 6.5 % most central events. 





4f) A ClpV/r 


80 A fipV/r 
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±0 03 % 

— 1 — KJ • KJ KJ / KJ 
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±0.02 % 


X 2 , vertex cut 
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-0.04 /0 


+0.13 (V 
-0.01 /o 
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-0.03 70 


SRC removal 


+0.07 0/ 
-0.04 /0 


+0.05 0/ 
-0.06 70 


+0.02<y 
-0.02 /0 


Total (no SRC removal) 


+0.42 co- 
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+0.18 (v 
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Total (with SRC removal) 
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-0.15 /o 
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-0.14 /o 



is not available for all data sets, we used TPC tracks only for all final results 
presented below and included the uncertainty into the systematic error, as- 
suming it is the same at all beam energies. 

The systematic contribution of track quality cuts was investigated by varia- 
tion of these cuts (x 2 , vertex) within reasonable limits. The observed effect on 
the reconstructed value of S PT was also included into the systematic error. 
We developped a strategy to remove from our data set contributions of short- 
range correlations (SRC) which may alter the measured fluctuation signal. 
The occurence of short-range correlations is due to quantum statistics, final 
state Coulomb interactions and the finite two-track resolution. The system- 
atic uncertainty introduced by this procedure is also listed in Table 2. For a 
detailed description of the procedure see Section 4. 

The sum of all systematic uncertainties reported in Table 2 is small compared 
to the experimentally observed values of S PT ~ 1%. 



4 Results on mean p T fluctuations 

We obtain a first qualitative hint for the presence of non-statistical mean p T 
fluctuations by an investigation of the event-by-event mean px distributions as 
shown in Fig. 3. The larger width at 40 A GeV/c is a consequence of the smaller 
azimuthal fiducial window and hence smaller multiplicity in this data set. We 
compare the event-by-event mean p T distributions to reference distributions 
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Fig. 3. Top : Event-by-event mean pr distributions in the 6.5% most central events 
at 40, 80, and 158 A GeV/c. Circles show the distributions of data events, solid lines 
indicate the mixed events. Bottom: Ratio between the distributions of data events 
and mixed events for 40, 80, and 158 A GeV/c. 



obtained by event mixing. The mixed events are constructed from particle 
momenta randomly chosen from data events of the same centrality class. Only 
one particle per measured event is used for a given mixed event, and the 
multiplicity distribution of mixed events is generated by sampling that of the 
data events. We calculated S PT and $ PT for the mixed event samples and 
found them to be consistent with zero within statistical errors at all three 
beam energies. 

The mixed event mean p T distributions exhibit a Gamma distribution shape [23] 
The subtle but clearly significant differences between the data and mixed event 
distributions are emphasized in Fig. 3 (bottom), where the ratio of the two 
is shown. The real event distributions are slightly wider, indicating a small 
but finite non-statistical contribution to the mean fluctuations at all three 
energies. A preliminary account of these results was presented in [24]. 
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Fig. 4. The two-track efficiency as function of the opening angle. 

Before we turn to a quantitative analysis of the observed non-statistical fluctu- 
ations, we discuss the effect of short-range correlations, which are not present 
in our mixed event sample but which may contribute to the fluctuation sig- 
nal. Physical short-range correlations are caused by Bose-Einstein quantum 
statistics and the final state Coulomb interaction; detector effects such as the 
finite two-track resolution lead to anti-correlations due to track merging. 

In Fig. 4 we present the opening angle distribution of particle pairs detected in 
the TPC divided by the corresponding mixed-event distribution. The resulting 
two-track detection efficiency is normalized to unity at large opening angles 
a. The slight increase towards smaller a can be attributed to small angle 
correlations such as Bose Einstein correlations and flow. At a < 10 mrad the 
efficiency drops due to track merging. 

It was demonstrated by NA49 [7] that the effects of Bose-Einstein correlations 
and two-track resolution on the fluctuation signal quantitatively cancel in their 
data set. In general, this is not necessarily the case but rather depends on the 
details of the momentum range under investigation, the bending power of the 
spectrometer and the two-track resolution. 

Short-range correlations show up at small momentum differences q and can 
be investigated by a study of the two-particle correlation function C^gmv)- 

The four-momentum difference q- mv = \J q 2 — q$ is the momentum difference 
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Fig. 5. The correlation function (^(qw) as function of q- mv for like-sign pairs (left 
panel) and unlike-sign pairs (right panel) at 158 A GeV/c. The histograms show 
raw data, the open circles are after removal of short-range correlations. The full 
dots show the reference correlation functions with a > 10 mrad (see text). 

in the pair rest frame, where q and q are the differences in three-momentum 
and energy of a particle pair assuming the pion mass for each particle. A 
systematic analysis of the two-pion correlation functions at SPS energies has 
been presented in [25]. 

The raw correlation functions C2,raw(<?inv) for pairs of like-sign and unlike-sign 
charged particles in 158 A GeV/c are shown as histograms in Fig. 5. The 
like-sign pairs show a positive correlation at q- mv < 70 MeV/c due to Bose- 
Einstein statistics in the pion-dominated sample. At very small qw, we observe 
an anti-correlation caused by the final state Coulomb repulsion and the finite 
two-track resolution. For the unlike-sign pairs, the positive correlation arises 
due to Coulomb attraction, and is again suppressed by the two-track resolution 
at small q- mv . As a reference, we also constructed the correlation functions 
C 2 (gi nv ) for signal and background pairs with opening angles a > 10 mrad 
(full dots in Fig. 5). These correlation functions are free of distortions from 
the finite two-track resolution but suffer from poor statistics at small q inv . 

Based on the experimental two-track resolution and the measured correla- 
tion functions C 2 (gi nv ), a procedure was applied to remove contributions from 
short-range correlations on the observed fluctuation signal. The procedure im- 
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plies a small modification of the measured events. To account for the two-track 
resolution, artificial tracks are added, and a small number of tracks from real 
pairs are removed to suppress the positive correlation at small q. After the 
procedure, the correlation functions C 2 (qmv) are flat, as indicated by the open 
symbols in Fig. 5. From the modified events the fluctuation signal is calcu- 
lated, which is then expected to be free of contributions from short-range 
correlations. The details of the procedure are described below. 

In a first step tracks are randomly added to form, with existing tracks, small 
opening angle pairs according to the observed two-track reconstruction effi- 
ciency The probability to add a track at a given a is determined from the 
distribution presented in Fig. 4 and constrained by the strength of the refer- 
ence correlation function of pairs with a > 10 mrad in Fig. 5. The fraction 
of artificially added tracks is a few percent, and their transverse momenta are 
randomly chosen from tracks of different events of the same centrality class. 

In a second step, positive correlations are eliminated by removal of tracks from 
small q inv pairs. The probability to remove a track is based on the measured 
reference correlation function C 2 ((/inv)- We define, for each pair of particles 
from the same event, the probability to reject randomly one of the particles: 

Prej = bll- n ] ) at q inv < 100 MeV/c, 

V ^2{qinv)J 



where C2(<2w) is the reference correlation function of pairs with a > 10 mrad 
and b is a parameter tuned for like-sign and unlike-sign pairs separately. After 
both steps, the obtained correlation functions are flat, as indicated by the 
open circles in Fig. 5, for typical values of b around 0.8-1.0. 

In Fig. 6 is presented the normalized dynamical fluctuation S PT in central 
Pb-Au events at 40, 80, and 158 A GeV/c as function of the size of the 
azimuthal acceptance window A0. The filled squares show the result for the 
unmodified events. If additional tracks are added to compensate for the two- 
track resolution, the results increase very slightly, as indicated by the open 
symbols. The filled dots show the results after the full SRC removal, which 
contains the compensation for the two-track resolution and the suppression of 
correlated small q pairs. The contributions from both effects do not cancel but 
lead to a reduction of S PT by 20 — 35% compared to the values obtained from 
unmodified events. The systematic errors of the correction procedure have 
been estimated to be less than 0.05% (absolute) in most cases by variation of 
the probabilities to add or reject a track within reasonable limits. They are 
small compared to the overall systematic uncertainties discussed before (see 
Table 2). 

Note that the azimuthal acceptance of the 40 A GeV/c data set is limited, 
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Fig. 6. The fluctuation measure S pT as function of the azimuthal acceptance window 
Acj) in central Pb-Au collisions at 40, 80, and 158 A GeV/c. The filled squares are 
without SRC removal, the filled circles are after the full SRC removal procedure. 
The open circles show the results obtained after compensation of the two-track 
resolution only. 



but no significant change of the fluctuation strength is observed when going 
from limited to full acceptance at the higher beam energies. This observation 
is the basis of the estimate of the systematic error introduced by the limited 
acceptance at 40 A GeV/c as discussed in Section 3.3. This implies also that 
azimuthal anisotropies of the particle distribution, such as elliptic flow, have 
little effect on the observed fluctuation strength. 

As a consistency check of the SRC removal procedure we applied the subevent 
method [19] to calculate the fluctuation S PT . If the subevents are based on a 
random selection of tracks from a given event, the observed fluctuation S PT 
is consistent with the result obtained by our standard method using Eq. (10) 
without SRC removal. The contribution from short range correlations can be 
estimated if the subevents are separated in pseudorapidity. We have chosen 
the pseudorapidity ranges 2.2 < rj < 2.4 and 2.5 < rj < 2.7 for subevents 1 
and 2, respectively. In this case, the observed mean px fluctuations are reduced 
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Fig. 7. The fluctuation measure S PT as function of the number of participants in 
Pb-Au collisions at 40, 80, and 158 A GeV/c. The data points are obtained after 
removal of short-range correlations. Also shown are extrapolations from pp data 
assuming N part scaling (dotted lines) and iV coll scaling (dashed lines, see text). 

compared to the case of random subevents, and consistent with those obtained 
using our previously described method for SRC removal. 



To facilitate a comparison of the data to models, which generally do not con- 
tain short-range correlations, only the results after SRC removal are shown 
below. A compilation of all results before and after SRC removal can be found 
in Table 3 and in Appendix A. 

In Fig. 7 the normalized dynamical fluctuation S PT is shown as function of the 
centrality of the collision at three different beam energies. The centrality is 
expressed in terms of the average number of participating nucleons (iV part ) in a 
given centrality bin. We observe that, at all energies, the dynamical fluctuation 
S PT is comparable in strength and decreases with centrality by approximately 
40% from about 1.3% to 0.7% over the centrality range under investigation. 

In this context, it is interesting to compare the present results to measure- 
ments in hadron-hadron collisions. In pp interactions particles are produced 
in a correlated way which leads to large non-statistical fluctuations. At the 
ISR, dynamical mean p t fluctuations have been measured in pp reactions at 
^/snn = 30.8 — 63.0 GeV [9]. Independent of beam energy, a value of 12% 
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was observed for S PT . In aa reactions, the observed dynamical fluctuation is 
reduced to about 9%. This was stated in [9] to be consistent with the larger 
average number of nucleon-nucleon interactions in aa of about 1.8. Under the 
assumption of an incoherent superposition of independent nucleon-nucleon 
collisions, S PT is expected to scale with (A^ co n): 

S P A T A = ■ <AUr 1/2 . (13) 



Use of this expression to extrapolate from the dynamical fluctuation of 12% 
observed in pp (N col \ = 1) to Pb-Au yields values of about 0.4% for the most 
central events. These extrapolated numbers are significantly below the data, 
as shown in Fig. 7, indicating that A-A collisions at SPS energies are not a 
straight superposition of nucleon-nucleon collisions also in this observable. 

A different approach is a possible scaling with charged particle multiplicity: 



Since the number of charged particles was found to scale close to linear with the 
number of participants (N paxt ) at SPS [26,27,28], we replace ((A / " pp )/(A / "aa)) 1/2 
by ((7V paJt )/2)-V2 : 

S A T A = m ■ «AW>/2)- 1/2 . (15) 



A comparison of Eq. (15) to our data is indicated by the dotted curve in Fig. 7. 
Although our data are at still somewhat smaller y / J/v r /v than explored in [9], 
we find that the magnitude of the fluctuation strength observed in Pb-Au is 
similar to the extrapolation from pp, if multiplicity scaling is assumed. 

In Fig. 8 we compare our mid-rapidity result for E PT in central Pb-Au collisions 
at 158 A GeV/c to the measurement from NA49 obtained at forward rapid- 
ity [7]. The NA49 result for $ PT [7] was used to evaluate S PT using Eq. (11). 
While the mean p T fluctuations measured at forward rapidity are consistent 
with zero, we observe a finite value of about 0.8% at mid-rapidity. Also shown 
are calculations from the RQMD [29] and URQMD [20] event generators. While 
the fluctuations observed in URQMD are small at all rapidities, RQMD exhibits 
a pronounced rapidity dependence which qualitatively describes the trend ob- 
served in the data. However, the RQMD absolute value around mid-rapidity is 
about 20% larger than that observed in the data. 

The magnitude of mean transverse momentum fluctuations from RQMD and 
URQMD was studied in more detail, as shown in Fig. 9. We compare our re- 
sult for S PT in central 158 A GeV/c Pb-Au collisions to calculations from 
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Fig. 8. The fluctuation measure S PT as function of pseudorapidity in central Pb-Au 
collisions at 158 A GeV/c. The CERES data are after SRC removal. Results from 
RQMD and URQMD calculations are shown as solid and dashed lines, respectively. The 
statistical errors of the calculations are indicated by the error bars at 77 = 2.45. The 
NA49 data point for the 5% most central Pb-Pb collisions (0.005 < p T < 1.5 GeV/c) 
is calculated from [7] (see text). The horizontal error bars indicate the rj acceptance 
range for the data points. Note that NA49 data were measured in the pion rapidity 
interval 4 < y w < 5.5. 

both models at 2.2 < rj < 2.7. Different mechanisms included in the mod- 
els have been investigated separately. If hadronic rescattering is switched off, 
both models show reasonable agreement with the data. In RQMD, the inclu- 
sion of rescattering leads to an increase of fluctuations. This behaviour is 
not understood and needs further investigation. In URQMD, the inclusion of 
meson-meson rescattering has very little effect on the observed fluctuations, 
however, inclusion of meson-baryon rescattering reduces the fluctuations con- 
siderably. If both rescattering modes are included in URQMD (labelled 'default' 
in the figure), URQMD falls below the data. We checked that also at 40 and 
80 A GeV/c the agreement between data and URQMD without rescattering 
is reasonable, while inclusion of rescattering in URQMD leads to a significant 
reduction of £ PT . At 40 A GeV/c, rescattering leads even to negative values 
for S PT . This indicates that secondary meson-baryon scattering has a strong 
effect on the observed fluctuation strength, and might be over-estimated in 
URQMD. It should be emphasized that the evaluation of event-by-event fluctu- 
ations may provide important information about the dynamics of the system 
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Fig. 9. Model calculations of T, PT from urqmd and RQMD in comparison to the 
CERES measurement (after SRC removal) at 2.2 < 7] < 2.7 in central Pb-Au 
collisions at 158 A GeV (see text). The results of the model calculations are shown 
for different rescattering scenarios. The dashed lines indicate the systematic error 
of the data. 

which can not be derived from inclusive observables. 



A compilation of the normalized dynamical fluctuation S PT measured at mid- 
rapidity and at different beam energies is shown in Fig. 10. The upper scale 
indicates the baryon chemical potential \xb at chemical freeze-out, related to 
\j $nn via a phenomenological parametrization given in [30]. The data point 
shown at ^/snn = 130 GeV is from the STAR experiment at RHIC [31] and 
measured over a somewhat wider pt range (0.1 < px < 2 GeV/c). A mea- 
surement by PHENIX [32] yields £ PT = 1.4l5J% and is consistent within its 
larger error. 2 At the SPS, the normalized dynamical fluctuation S PT is about 
0.7%, while larger fluctuations of 1.2% are observed at RHIC. The STAR data 
are not corrected for the contribution of short-range correlations, however, the 
effect was estimated to be small (less than 10%) [31]. The evolution of S PT 
with beam energy looks smooth and does not show any indication of unusually 
large fluctations at any beam energy. 

Models predict enhanced mean p T fluctuations if the system has passed close to 

2 The measure Sp T was calculated using the PHENIX measurement for d and uip T 
from [32] and the approximation £ 2 ~ 2dw PT . 
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Fig. 10. The fluctuation measure Xp T as function of ^^saT/v and of \ib &t chemical 
freeze-out [30]. The full circles show CERES results (after SRC removal) in central 
events at 40, 80, and 158 A GeV/c. The brackets indicate the systematic errors. 
Also shown is the STAR result [31] at y / ijv r /v = 130 GeV which is not corrected 
for SRC. Results and statistical errors from rqmd and URQMD calculations (with 
rescattering) are indicated as solid and dashed lines, respectively. 

the critical point of the QCD phase diagram. At SPS energies and for the finite 
rapidity acceptance window of the CERES experiment, the fluctuations should 
reach values of about 2%, i.e. more than three times larger than observed in the 
present data 3 . Most important, no indication for a non-monotonic behaviour 
as function of the beam energy has been observed. This suggests that the 
critical point may not be located in the \xb regime below 450 MeV. 

The results from RQMD and URQMD show rough agreement with the data, 
except for the URQMD calculation at 40 A GeV/c where S PT is negative (see 
Fig. 10). We note that a positive value of S PT = 0.38^o;4g% is obtained from 



3 The predicted fluctuations in the measure 
2% in S PT in the CERES acceptance [33]. 



F = 1.1 in [13] corresponds to about 



22 



Beam momentum (A GeV/c) 


/I n 
4U 


on 

80 


loo 


S PT (%) 
(no SRC removal) 
(with SRC removal) 


i.uo ± u.uy_Q 

n 71 +0.12+0.42 
u - ' i -0.15-0.14 


n Q/1 +0.08+0.18 
^^-0.09-0.14 

fl Rc+0.11+0.19 
U - DO -0.14-0.15 


n qq+0. 06+0.15 
u - yy -0. 07-0.14 

U. to ± U.Uo_q 14 


<& PT (MeV/c) 
(no SRC removal) 


Z.4 ± U.4_j 2 


Q O 4- n 7+1-0 

3.9±0.7_ 06 


O.U ± U.D_ 1 2 


l^witn oxiL/ removal i 


ii i n a +1-5 
1.1 ± 0.4_ x 2 


i o i n 'r+l.l 

l.y ± U. i _q n 


q 1 4- n ^+2. 1 
o.l ± U.D_^ 2 


n 




loDZ I Z 






4 7 _i_ n m 


± ,\J ± _1_ U.Ut: 
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438.67 ± 0.05 


434.68 ± 0.06 


425.20 ±0.05 


^Ap^ (MeV/c) 


290.14 ± 0.05 


284.17 ±0.06 


278.97 ±0.05 


^/(AM^) (MeV/c) 


37.919 ±0.006 


25.642 ±0.005 


22.613 ±0.004 



Table 3 

Summary of mean pt fluctuations in the 6.5 % most central Pb-Au events at 40, 
80, and 158 A GeV/c. Statistical and systematic errors are quoted for S PT and & PT . 
The other quantities are without SRC removal and statistical errors only are shown. 



URQMD at 40 A GeV/c, if rescattering is switched off. 

Mean transverse momentum fluctuations can be related to event-by-event fluc- 
tuations of the temperature. The following expression was given in [34]: 

2 

% T = V2 (N) (16) 



where (T) is the average temperature which fluctuates event by event with a 
standard deviation u T . To derive an upper limit for the temperature fluctu- 
ations, it is assumed that the full magnitude of the observed mean px fluc- 
tuations after removal of short-range correlations is due to temperature fluc- 
tuations. For a temperature T = 160 MeV, values for ot of 1.4±0.3 MeV, 
1.3±0.3 MeV, and 1.5±0.2 MeV are obtained in the 6.5% most central events 
at 40, 80, and 158 A GeV/c, respectively. For T = 120 MeV, the results are 
smaller by about 15%. We conclude that event-by-event temperature fluctua- 
tions are at most 1% of the average temperature. Such small numbers point to 
a large degree of thermalization of the system. We note that small temperature 
fluctuations are expected also in the case of complete thermalization [11,34] 
and may in that case be used to determine the heat capacity of the system [34]. 
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5 Conclusions 



We have presented event-by-event fluctuations of the mean transverse mo- 
mentum in Pb-Au collisions at 40, 80, and 158 A GeV/c. In central events 
at these beam energies, significant non-statistical positive mean p T fluctua- 
tions are observed. Based on the measured two-particle correlation functions 
and the experimental two-track resolution a procedure was derived which ac- 
counts for the contributions of short-range correlations present in the data. 
This procedure reduces the observed mean pr fluctuations by 20-35% but they 
remain significantly different from zero. At all beam energies, we find for the 
normalized dynamical fluctuation S PT values of about 0.7% of the mean p T . 
These results at mid-rapidity are complementary to previous measurements at 
158 A GeV/c in the forward hemisphere (4 < y n < 5.5) where no significant 
non-statistical fluctuations were found [7]. 



The present mid-rapidity results at SPS are somewhat smaller than the val- 
ues measured by STAR at y/s NN = 130 GeV. The existing data from different 
beam energies do not exhibit an indication for unusually large fluctuations or 
a non-monotonic behaviour, which might have pointed to the crossing of the 
critical point of the QCD phase diagram somewhere in this range of beam en- 
ergy or baryon chemical potential [13]. Further studies at lower beam energies 
are needed for a full exploration of the QCD phase diagram. 



At all three beam energies under investigation we observe a systematic de- 
crease of the mean pt fluctuation strength with increasing centrality of the 
collision. Based on measurements of non-statistical mean px fluctuations in 
pp-collisions at the ISR we find that the centrality dependence of mean p? 
fluctuations in Pb-Au is consistent with an extrapolation from pp, assuming 
that mean p T fluctuations scale with (iV) -1 / 2 . This may suggest that the fluc- 
tuation pattern is not reduced strongly by potential rescattering of hadrons 
after hadronization. A comparison to results from the RQMD and URQMD mod- 
els indicates that secondary rescattering, if enabled, tends to decrease the 
fluctuation strength, while calculations without rescattering show reasonable 
agreement with the data. Such a scenario is supported by the observation 
of high densities and a short mean free path at thermal freeze-out, as de- 
rived from a recent analysis of pion interferometry data [35], which point to 
a short lifetime of the hadronic phase. This issue needs further investigations 
of smaller collision systems and minimum bias Pb-Au data, as well as more 
detailed model studies. 
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A Summary of mean p T fluctuations 



In this appendix the results of mean p T fluctuations in Pb-Au collisions at 
40, 80, and 158 A GeV/c and different centrality classes are summarized. The 
results for S PT and $ PT are listed before and after removal of short range 
correlations (SRC). They contain statistical and systematic errors. All other 
parameters are before SRC removal and statistical errors only are shown. 

Table A.l 

Summary of mean pt fluctuations at 40 A GeV/c. 

Centrality - 5 % 5 - 10 % 10 - 15 % 15 - 20 % 

S PT (%) 

(r\n SRC! rpmoval^ i n4 +0.10+0.51 1 -.0+0.12+0.33 ,.,+0.12+0.31 , fi o+0.16+0.32 

(no snu removal; i-U4_ n _ 26 J-J-8_0.13-0.22 ^^-o.is-o.ie i - D8 -o.i8-o.25 

fwitli STIP rpmnvan D «Q+0.14+0.52 -n+0.17+0.34 -, o-, +0.15+0.31 -, 9fi +0.21+0.32 

(witn situ removal; u.oy_ 017 _ 027 u.^_ 21 _ 23 i._i_ 17 _ 16 i-^o_ 26 _ 025 
% T (MeV/c) 

(no SRC removal) 2.3 ± 0.4±g;| 2.4 ± 0.5±J;| 3.3 ± 0.5±J;t 3.3 ± 0.6±J;| 

(with SRC removal) 1.1 ± 0.4±g;| 1.1 ± 0.5±J;| 2.0 ± 0.5±J;| 1.9 ± 0.6±J;| 

n 495951 396868 327473 214202 

(N) 61.19 ±0.01 49.39 ±0.01 40.94 ± 0.01 33.88 ± 0.01 

pT(MeV/c) 439.03 ±0.05 436.07 ± 0.07 433.15 ±0.08 429.75 ±0.11 



Ap^ (MeV/c) 290.35 ±0.05 288.71 ± 0.07 286.90 ± 0.08 284.84 ±0.11 



(AM 2 ) (MeV/c) 37.395 ± 0.007 41.400 ± 0.009 45.333 ± 0.012 49.464 ± 0.018 
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Table A.2 

Summary of mean px fluctuations at 80 A GeV/c. 



Centrality - 5 % 5 - 10 % 10 - 15 % 15 - 20 % 



fnn qRP rpmnvan n qi+ll-^+U-^ 1 n4 +U.uy+U.l I -, 97 +U.lU+U.l c -, c-q+U.lb+U.iy 

(no brt^ removal J U.y4_ 010 _ 026 i-U4_ 010 _ 018 i-^ ' _ .n-o.i4 i - Dd -o.i8-o.i9 

19 
19 



(with SRC removal) 0.70±°;^J 8 7 0-78^^ 0.99±°;^ l^l^ 19 



<S> PT (MeV/c) 

(no SRC removal) 4.0 ± 0.St\i 3.9 ± 0.7l?;{j 5.0 ± 0.8±g;^ 6.5 ± 1.4+^ 

(with SRC removal) 2.3 ± 0.8^{^ 2.2 ± 0.7l?;{j 3.1 ± 0.8±J;2 4.6 ± 1.4^ 

n 135075 161036 137261 45647 

(N) 129.22 ±0.05 110.20 ±0.04 92.22 ± 0.04 82.19 ±0.07 

pf (MeV/c) 434.95 ± 0.07 433.12 ± 0.07 430.55 ± 0.08 428.90 ± 0.15 



yX^4 (MeV/c) 284.26 ± 0.07 283.57 ± 0.07 282.15 ± 0.08 281.27 ± 0.15 
yj (AMp T ) (MeV/c) 25.340 ±0.006 27.384 ±0.007 29.888±0.008 31.710±0.016 



Table A.3 

Summary of mean pt fluctuations at 158 A GeV/c. 



Centrality - 5 % 5 - 10 % 10 - 15 % 15 - 20 % 

S PT (%) 

(no SRC removal) 0.98l°;°£t°;!| l.lllg^llue L38 ± 0.08±g;^ 1.45 ± 0.10±g;l? 

(with SRC removal) C^gggg 0.90±g;gg±g;^ 1.12±g;gg±git 1.16±g;jg±g;^ 
$ PT (MeV/c) 

(no SRC removal) 5.1 ±0.7+^ 5.4±0.7l?;f 6.8 ± 0.8±?;| 6.5 ± 0.9±?;| 

(with SRC removal) 3.3 ± 0.7±} | 3.6 ± 0.7l?;f 4.4 ± 0.8l?J 4.1 ± 0.9±?;| 

n 151191 153713 137414 102787 

(AO 161.85 ±0.05 135.35 ±0.05 113.06 ± 0.05 97.47 ± 0.05 

pf (MeV/c) 425.40 ± 0.06 423.76 ± 0.06 421.70 ± 0.07 419.49 ± 0.09 

^/Apf(MeV/c) 279.06 ± 0.06 278.30 ± 0.06 267.23 ± 0.07 276.18 ±0.09 

,/(AM p 2 T ) (MeV/c) 22.335 ±0.005 24.379 ±0.005 26.723±0.007 28.612±0.009 



28 



